%This function calculates the junction area as a function of the width of
%the tape, w, and the junction angle, phi
function [Areaslist]=areacalculation(w,testangles)
Areaslist=zeros(1,length(testangles));
for i=1:length(testangles)
    %w=12.7E-2; %mm
    testangle=testangles(i);
    Area=@(phi) (1./2).*(w.^2).*((secd(phi./2)).^2).*(cotd(phi./2));
    Areaslist(i)=Area(testangle);%ones(1,length(testangle));%
end